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INTRODUCTION 


The  long  range  objective  of  the  research  presented  in  this  report 
was  to  deVelpp  a  basic  understanding  cf  low  velocity  detonation  (LVD)  and 
to  determine  conditions  for  its  initiation  and  propagation.  Such  knowledge 
is  required  for  predicting  the  LVD  behavior  of  energetic  materials  used 
in  propulsion  systems.  It  is  important  to  the  Air  Force  because  the 
possibility  of  initiating  LVD  with  pressures  below  10  kbars  constitutives 
a  major  hazard  in  the  efficient  use  and  manufacture  of  liquid  jnonopropel- 
lants  and  energetic  plasticizer  compounds.  Events  that  must  be  safeguarded 
against  are  the  malfunction  of  rocket  motors  caused  by  LVD  in  the  fuel 
container  or  fuel  feed  system,  and  serious  accidents  caused  by  LVD  during 
the  preparation  and  transportation  of  large  quantities  of  such  materials. 

1—6 

Shock  wave  studies  .during  the  last  decade,  extending  the  earlier 

7 

work  discussed  by.  Taylor,  have  established  that  most  liquid  propellants 
and  explosives  exhibit  two  mGdes  of  self-sustaining  detonation.  These 
modes  have  distinct  velocities  of  propagation  and  have  been  characterized 
accordingly  as  high-velocity  detonation  (HVD)  and  low-velocity  detonation 
(LVD).  HVD  propagates  at  approximately  6  mm//Usec  with  a  pressure  in  the 
100-kilobar  region  and  can  be  initiated  by  a  shock  with  a  pevak  pressure 
of  about  60  kilobars.  LVD  propagates  with  a  slightly  higher  velocity 
than  ambient  sound  speed,  at  approximately  2  mm/ /Usee  with  a  pressure  in 
the  10-kilobar  region,,  und  can  be  initiated  by  a  shock  with  a  peak 
pressure  of  about  4  kilobars.  in  some  cases,  however,.  LVD  has  been 
observed  to  transform  into  HVD. 
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Since  explosives  are  required  to  behave  reproduclbly  and  produce 
pressures  in  the  100-kilobar  region  in  most  practical  applications, 
detonation  studies  have  concentrated  on  HVD  rather  than  on  LVD.  Con¬ 
sequently,  conditions  for  initiation  and  propagation  of  LVD  are  not  as 
well  established  as  those  for  HVD,  and  our  understanding  of  LVD  is  not 
as  complete  as  'that  of  hVl*.  This  situation  is  further  accentuated  by 
the  disparity  in  the  theoretical  treatments  of  HVD  and  LVD.  The  hydro- 
theriuddyhamic  theory  of  detonation  Is  satisfactory  for  modeling  the 
initiation  and  propagation  of  HVD  and  for  .caleulaiing-JIVD  detonation 
parameters  with)  the  Chapman- Jouguet  (C-J)  hypothesis  when  the  equation 
of  state  of  the  i detonation,  products  is  known,  but  there  is  no  comparable 
theory  for  LVD.  ... 


BACKGROUND  AND  APPROACH 


Previous  research  on  LVD  in  Poulter  Laboratory  at  Stanford  Research 

Institute  (SRI)  was  supported  by  AFOSR  under  Contract  F44620-69-C-CQ79  and 
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is  documented  in  the  final  report  written  for  this  contract.  This 

research  included  both  theoretical  and  experimental  studies.  Computer 

calculations  were  performed  in  the  theoretical  study  to  test  the  cavita- 
5 

tion  model  of  LVD,  and  ytterbium  stress  gages  were  used  in  the  experi¬ 
mental  study  to  determine  the  states  attained  in  LVD  waves.  A  brief 
discussion  of  this  work  is  given  below  as  background  for  the  present 
'project. 

In  the  cavitation  model,  LVD  is  treated  as  a  reactive  shock  propa¬ 
gating  in  a  cavitated  liquid.  The  liquid  is  first  cavitated  by  inter¬ 
actions  with  the  container  wall,  and  reaction  is  then  initiated  by 
bubble  collapse  in  the  compressed  cavitation  field.  Computer  studies 
were  performed  to  investigate  the  cavitation  process  and  to  determine 
the  amount  of  reaction  required  to  support  LVD.  Two-dimensional  computer 
calculations  of  a  shock  propagating,  in  a  liquid  contained  in  a  cylindrical 
steel  tube  showed  that  precursor  wall  waves  produce  regipns  of  tension 
sufficient  to  cavitate  the  liquid  before  it  is  compressed  by  the  main 
shpck.  Simulation  of  hot-spot-initiated  reaction  in  the  calculations  by 
addition  of  energy  behind  the  shock  leu  to  the  conclusion  that  only  part 
of  the  chemical  energy  is  required  to  support  the  wave.  Methods  of  calcu¬ 
lating  LVD  parameters  were  formulated  with  the  hydrothermodynamics  of 
the  cavitation  model,  and  a  complete  equation  of  state  for  nitromethane 
was  introduced  into  the  TIGER  code  so  that  detonation  calculations  could 
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be  made  for  incomplete  reaction.  The  calculations  showed  that  only  about 
20  percent  of  the  liquid  is  required  to  react  to  support  LVD  and  that  LVD 
parameters  are  insensitive  to  the  initial  degree  of  cavitation.  The 
results  of  this  theoretical  study  are  significant  because  they  place  the 
cavitation  model  on  a  firm  foundation  and  justify  its  use  as  a  basis  for 
further  studies  to  J^velop  a  satisfactory  treatment  o'*  LVD. 

Ytterbium  gages  were  used  to  study  LVD  In  ethyl  nitrate  and  in  FEFO. 

LVD  propagation  velocities  were  found  to  be  about  2  mm/jisec  in  ethyl 
nitrate  arid  3  mm/jUsee-  in  FEFO;  peak  pressures  were  recorded  from  4  to  7 
kbars. 

The  present  work  on  LVD  was  a  logical  continuation  of  that  discussed 
above.  Theoretical  and  experimental  studies  were  continued  to  determine 
conditions  for  the  onset  and  termination  of  reaction  in  the  compressed 
cavitation  field  and  conditions  for  thei  initiation  of  LVD  and  its  transition 
into  HVD. 

Two  approaches  were  taken  in  the  theoretical  study.  The  first  approach 
was  8.  numerical  study  of  bubble  dynamics  with  heat  and  mass  transfer  and 
chemical  reaction.  The  idea  was  to  try  to  establish  conditions  for  the 
onset  and  termination  of  reaction  in  a  single  bubble,  >and  then  to  try  to 
apply  these  conditions  to  the  cavitation  field  to  determine  the  amount 
of  reaction  supporting  LVD.  The  computer  calculations  on  a  single  bubble 
were  performed  by  J.  Levine  and  D.  C.  Wooten  at  Ultrasystems,  Inc.,  and 
an  account  of  this  work  is  given  in  the  Appendix.  The  second  approach, 

i 

more  empirical  than  the  first,  was  to  study  reactive  shock  discontinuities. 
Reactive  Hugrmiot  curves  were  constructed  in  the  pressure-volume  (p-v) 
plane  to  account  for  the  initiation  and  propagation  of  LVD  and  for  its 
transition  into  HVD.  These  approaches  will  lead  to  a  satisfactory  model 
for  LVD  when  the  parameters  used  in  the  treatment  of  reactive  shock 
discontinuities  can  be  evaluated  with  the  results  of  the  computer 
calculations,  and  these  parameters  agree  with  experimental  values. 
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A  new  method  for  placing  multiple  stress  gages  in  a  liquid  was 
devised  in  the  experimental  study.  This  new  method  of  gage  emplacement 
allows  the  gages  to  move  with  the  flow  and  record  a  more  accurate  and 
quantitative  description  of  the  states  attained'  in  LVD.  Experiments 
using  multiple  ytterbium  stress  gages  were  performed  on  mixtures  of 
nitromethane  and  tetranitromethane. 


THEORETICAL  STUDY  OF  LVD 


LVD  Treated  as  a  Reacting  Shock  Discontinuity 

Although  a  complete  treatment  of  LVD  must  account  for  liquid- 
container  wall  interactions  and  the  dynamics  and  thermodynamics  of  a 
cavitation  field  with  heat  transfer,  mass  transfer,  and  chemical  reaction 
the  present  work  is  based  on.  the  assumption  that  a  detailed  account  of 
these  processes  is  not  required  to  model  the  initiation  and  propagation 
of  LVD  and  its  transition  into  HVD.  This  approach  is  motivated  by  the 
one-dimensional  CJ  model  for  gaseous  detonation,  which  is  satisfactory 
for  calculating  the  average  velocity  of  HVD  without  treating  the  complex 
system  of  interacting  transverse  waves  necessary  for  its  propagation. 


Detonation  will  be  considered  as  a  reactive  shock  discontinuity  and 
the  flow  behind  it  as  one  dimensional.  The  time  scale  for  compression 
and'  reaction  supporting  the  shock  is  therefore  assumed  to  be  very  small 
compared  with  the  time  scale  for  the  subsequent  pressure  release  in  the 
wave.  Consequently, ' HVD  is  represented  as  a  shock  discontinuity  with 
complete  reaction,  and  LVD  as  a  shock  discontinuity  with  incomplete 
reaction.  Let  P,  v,  p,  and  e  denote  density,  specific  volume, 
pressure,  and  specific  internal  energy,  and  let  the  subscript^  H  and 
o  denote  quantities  immediately  behind  and  in  front  of  a  6hock  dis¬ 
continuity.  Then  the  Rankine-Hugoniot  jump  conditions  expressing  the 
conservation  of  mass,:  momentum,  and  energy  across  a  shock  discontinuity 
propagating  into  stationary  material  can  be  written  as 


P  „v 
H  H 


2(e, 


H 


"  (”-VUH  “  VhVVc 

(1) 

-•«>  ■  Vvo-V  *  *5 

(2) 

i  Preceding  page  blank 
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where  U  denotes  the  shock  velocity,  u  denotes  the  particle  velocity, 

and  the  initial  pressure  p  has  been  neglected.  Let  t  and  r  denote 

o 

time  and  distance;  then  the  equations  expressing  the  conservation  of  mass, 
momentum,  and  energy  in  the  flow  behind  the  shock  can  be  written  as 

dp  Bu 

Tt  +  p3?  *  0  <3> 

P£  ♦  |£  =  0  (4) 

dt  or 


de  dv 

dt  y  dt 


where  d/dt  =  3/Bt  +  u3/3r  and  the  subscripts  on  the  partial  derivatives 
have  been  omitted  for  notational  convenience.  When  the  flow  is  supported 
by  an  exothermic  reaction  described  by  a  reaction  coordinate  X,  combina¬ 
tion  of  Eq.  (5)  and  the  e  =  e(p,  v,  X)  equation  of  state  leads  to  the 
equation 

dp  a  dP  IV  dX 

—  =  c  —  +  —  Q  —  (6) 

dt  dt  v  dt 

where  c  =  (3p/3p)^  denotes  the  sound  speed,  s  denotes  the  specific 
s  f  A 

entropy,  T?/v  -  (Bp/de)  .  is  the  Grlineisen  parameter,  and  Q  =  -(de/dX) 

v ,  A  P; 

denotes  the  specific  heat  of  reaction. 

To  eliminate  equation-of-state  problems  it  is  assumed  that  the 
explosive  mixture  behaves  as  a  polytropic  medium  described  by  an  e(p,v,X) 
equation  of  state  of  the  form 

e  =  e°  -  XQ  +  7-“-  (7) 

x  K  -  1 


where  K  denotes  the  polytropic  index  and  e denotes  the  heat  of 
formation  of  the  explosive.  Combination  of  Eqs.  (1)  and  (2)  with  Eq.  (7) 
leads  to  the  following  equations  for  a  strong  reactive  shock  discontinuity 
in  a  poly tropic  medium, 
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2  Uu  =  (K  +  l)Ua  +  2(K  -  1)  QX 
H  H  H 

2  p iv  =  (K  +  l)u *  +  2(K  -  1)  QX 

Ho  n  H 


(S) 

(9) 


Let  Up  denote  the  detonation  velocity  of  the  CJ  discontinuity, 

and  let  the  -subscript  CJ  denote  the  Chapman-Jouguet  state.  Then, 

substituting  X  -  1  and  the  CJ  condition  written  as  U  =  (K  +  i)u 

CJ  D  CJ 

into  Eqs.  (8)  and  (9)  yields  the  following  well-known  expressions  for 
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HVD, 

U2  =  2<K2  -  1)Q  (10) 

D 


p  =  2(K  -  1)QP 

CJ  o 


(11) 


The  equations  for  a  nonreactive  shock  discontinuity  are  obtained  by 
setting,  X^  =  0  in  Eqs.  (8)  and  (9).  It  follows  from  Eqs.  (8)  and  (9) 
that  the  release  of  chemical  energy  in  a  shock  discontinuity  lowers  the 
values  of  shock  pressure  and  particle  velocity  with  respect  to  the  values 
across  a  nonreactive  discontinuity  propagating  at  the  same  velocity. 

In  other  words,  the  reactive  Hugoniot  curve  lies  above  the  nonreactive 
Hugoniot  in  the  (p-u)  plane. 


The  differential  equation  governing  the  propagation  of  a  shock  dis¬ 
continuity  will  now  be  derived.  The  assumption  will  also  be  made  that 

the  flow  behind  the  shock  is  nonreactive.  Let  D(  )  /Dt  =  d(  )/dt  + 

H  H 

Ud(  )/dr  denote  the  directional  time  derivative  along  the  shock  path; 
H 

then  the  shock  pressure  and  particle  velocity  satisfy  the  following 
identities 


dp 

5p 

Dt 

=  + 
dtH 

<U  -  UH) 

Sr 

H 

% 

du 

(U  -  u  ) 

Su 

=  —  + 

SrH 

Dt 

dtH 

H 

(12) 


(13) 
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The  combination  of  Eqs.  (12)  and  (13)  with  Eqs.  (3)  and  (4),  and  Eq.  (6) 

subject  to  the  condition  dX/dt  =  0,  gives  the  required  equation  along 

H 

the  shock  path  as 


Dp  Du 

H  .  „  H 

- +  P  U - 

Dt  o  ot 


u  )a 
H 


a\  3u 
°H/  3r„ 


(14) 


At  this  stage  the  subscript  H  will  be  omitted  for  notational  simplicity 
because  we  will  be  dealing  only  with  quantities  evaluated  at  the  top  of 
the  shock  discontinuity. 

Our  consideration  of  reactive  shock  propagation,  based  on  Eq.  (14), 

is  an  extension  of  some  ideas  developed  for  Ballistic  Research  Labora- 

10 

tories  under  Contract  No.  DA-04-200-AMC-2469(X) .  The  particle  velocity 
gradient  at  the  shock  front  is  always  assumed  to  be  positive  because  this 
condition  is  satisfied  in  most  initiation  experiments.  Reactive  dis¬ 
continuities  satisfy  the  conditions  X  ^  0  if  U  <  U^,  and  X  =  1  if 
U  DX/dt  >  0  for  accelerating  shocks,  DX/Dt  <  0  for  decelerating 

shocks,  and  DX/dt  =  0  for  shock  propagating  at  constant  velocity.  It 
is  evident  on  physical  grounds  that  shocks  satisfy  the  condition  Dp/Du  >  0. 
Consequently,  Dp/Dt ■>  0  and  Du/Dt  >  0  for  accelerating  shocks, 

Dp/Dt  <  0  and  Du/Dt  <  6  for  decelerating  shocks,  Dp/Dt  =  Du/Dt  =  0 
for  constant  velocity  shocks,  and  reactive  Hugoniot  curves  satisfy  the 
conditions  Dp/DX  >  0  and  Du/Da  >  0  since  p  and  u  are  functions 
of  X.  Combination  of  these  conditions  with  Eq.  (14)  leads  to  the 
conclusions  that  (('.  ••  u)8  -  c2  >  0  for  accelerating,  shocks,  (U  -  u)2  - 
c  <0  for  decaying  shocks,  and  that  (U  -  u)  -  c  =0  for  constant 
velocity  shocks.  With  the  notation  introduced  early,  for  example,  the 

sonic  condition  expressed  by  the  last  equation  becomes  U  =  u  +  c 

D  CJ  CJ 

for  HVD,  jt  is  also  convenient  to  denote  the  propagation  of  LVD  by  U*. 
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Conditions  at  a  sonic  point  are  readily  obtained  by  rewriting  the 


term  in  brackets  on  the  right-hand  side  of  Eq.  (14)  as 


C(U  -  u)2  -  c2] 


i\/ucj\2 

lAVV  • 


(15) 


Equation  (15)  gives  X(u  /u)2  =1  as  the  condition  for  a  constant 

CJ 

3  2 

velocity  reactive  shock  discontinuity.  Since  (U  -  u)  -  c  =0  at  a 

point  where  a  Rayleigh  line  is  tangent  to  a  Hugoniot  curve.  >(»»  /u)2  =  1 

CJ 

is  also  the  condition  for  a  Rayleigh  line  to  be  taKf^nt  to  a  frozen 

Hugoniot  curve  in  the  (p-u)  plane.  In  other  words,  X(u  /u)  =  1 

CJ 

defines  the  locus  of  sonic  states  in  the  (p-u)  plane.  The  locus  of 

sonic  states  in  the  (p-v)  plane  is  readily  found  to  be  v  =  [Ky  /(K  +  l], 

o 

It  follows  from  Jouguet's  Ruie11  that  accelerating  shocks  satisfy 

the  condition  v  >  [Ky  /(K  +  1]  and  that  decelerating  shocks  satisfy 
H  o 

the  reverse  inequality  v  <  (Ky  /(K  +  l).  Consequently,  in  the  (p-v) 

H  o 

plane,  accelerating  shocks  are  represented  by  points  that  lie  to  the 
right  of  the  sonic  locus,  decelerating  shocks  are  represented  by  points 
that  lie  to  the  left  of  the  sonic  locus,  and  constant  velocity  shocks  are 
represented  by  points  lying  on  the  sonic  locus. 


We  are  now  in  a  position  to  consider  initiation  of  detonation  when 


the  initial  particle  velocity  u  <  U  ,  and  the  initial  pressure 

i  D 

p  <  p  .  Since  weak  shocks  do  not  promote  detonation,  it  is  assumed 
i  CJ 


that  shocks  are  unreactive  until  a  critical  pressure  p 

ip* 
l  ' 

* 


So  X  =  0  for  Pi  t 


i 

and  shocks  with  initial  pressures  p  i  p 


is  attained. 

* 


lie  on  the  unreactive  Hugoniot  curve.  Let  v 


i  *i 
denote  the  volume  on 


the  unreactive  Hugoniot  when  p^  =  p^  ,  and  consider  the  locus  of  initial 


states  when  p^,  is  increased  from  p^ 


to  p 


CJ 


As  p^  approaches 


p  ,  the  initial  states  must  lie  to  the  right  of  the  sonic  locus  .in  the 
CJ 

(p-v)  plane  because  shocks  with  p  <  p  promote  detonation.  Moreover 

i  CJ 


ii 


the  initial  states  nust  lie  to  the  left  of 


* 

as  approaches  p^  , 

jc  * 

the  sonic  locus  because  (p^  ,v^  )  lies  on  the  unreacted  Hugoniot  curve. 

It  follows  from  these  observations  that  the  initial  states  lie  on  a  curve 

connecting  (p  *,v  *)  tc  (p  ,v  )  and  that  this  curve  must  intersect 
i  i  CJ  CJ 

the  locus  of  sonic  points.  Let  p*  denote  the  pressure  at  this  point 
of  intersection  and  let  X*  denote  the  corresponding  value  of  the 
reaction  coordinate.  Then  reactive  shocks  with  initial  pressures  in 
the  range  p^*  <  p*  will  decay,  and  shocks  with  initial  pressures  in 
the  range  p*  <  p  <  p  will  build  up  to  HVD.  Furthermore,  it  is 

X  UJ 

possible  that  shocks  initiated  with  p^  =  p*  will  propagate  at  con¬ 
stant  velocity  and  that  the  sonic  point  will  represent  LVD.  But  the 
stability  of  such  shocks  must  be  considered  before  it  can  be  concluded 
that  this  is  indeed  the  case. 

10 

Conditions  for  the  initiation  of  HVD  with  p  =  p  lead  to  the 

i  CJ 

conclusion  that  a  sonic  point  will  be  stable  and  represent  LVD  when 
initial  shock  pressures  greater  than  this,  sonic  pressure  do  not  promote 
further  reaction.  Shocks  with  initial  pressures  higher  than  the  sonic 
pressure  are  therefore  represented  by  points  on  the  frozen  Hugoniot  curve 
passing  through  the  sonic  point.  Shocks  with  pressures  lower  than  the 
sonic  pressure  are  represented  by  points  on  a  reactive  Hugoniot  curve. 

Consider  the  case  when  there  is  only  one  sonic  point,  say 

(Pi,  CKv  /(K+l)],  and  let  X*  be  the  corresponding  value  of  the 
o 

reaction  coordinate.  Then  the  Hugoniot  curve  exhibits  two  branches: 

the  lower  branch  with  0  £  X  s  X*  and  p*  £  p_^  £  p*  has  a  positive 

slope;  the  upper  frozen  branch  with  X  =  X*  and  p^  ^  p*  has  a  negative 

slope.  LVD  is  represented  by  the  sonic  point  (p*,  X*).  Consequently, 

%  * 

shocks  with  p  in  the  range  p  £  p  <  pi  decay,  shocks  initiated 
1  11 

with  p^>  pf  decay  to  LVD,  and  those  initiated  with  p^  =  p*  propagate 
unchanged  as  LVD  with  a  velocity  U^. 


12 


Now  consider  the  case  when  the  Hugoniot  curve  has  a  similar  shape, 
a  positively  sloped  lower  branch  and  a  negatively  sloped  upper  branch, 
but  intersects  the  sonic  locus  in  two  points  (p*,  X*)  and  (pi,  X2) 
with  pz  >  p*  and  X2  >  X*.  In  this  case  the  upper  sonic  point  repre¬ 
sents  LVD.  Shocks  with  p,  in  the  range  p*  ^  p^,  <  p*  decay,  shocks 

with  p^  >  pa  decay  to  LVD,  shocks  with  p^  in  the  range  p*  £  p^  <  pj 

build  up  to  LVD,  and  shocks  with  p^  =  pa  propagate  at  constant  velocity 

as  LVD. 

As  indicated  by  the  computer  studies,  of  bubbles,  such  a  situation 

would  occur  in  the  cavitation  model  when  the  reaction  initiated  by 

bubble  collapse  is  terminated  by  bubble  expansion  and  does  not  go  to 

completion.  Relationships  for  Hugoniots  exhibiting  these  two  branches 

were  formulated  because  a  Hugoniot  of  this  type  was  observed  by  Evans 
12 

et  al.  in  a  shock  initiation  study  of  low  density  pressings  of  ammonium 
perchlorate  sponsored  by  AF03R  under  Contract  AF49(638)-1124. 

The  transition  of  LVD  into  HVD  is  accounted  for  by  assuming  that 

there  is  another  critical  pressure  p^y  for  the  onset  of  bulk  reaction. 

/ 

In  this  case  reaction  initiated  by  bubble  collapse  at  the  pressure  p^ 

is  not  terminated  at  later  times.  Bulk  reaction,therefore  occurs  at 
/ 

the  pressure  p^  .on  the  frozen  Hugoniot  curve,  .passing  through  the  sonic 

point  representing  LVD.  The  locus  of  shocked  states  with  p  in  the 

1 

range  p  ^  p  ^  p  forms  a  curve  that  again  intersects  the  sonic 
i  i  CJ 

locus  and  has  a  positive  slope  in  the  neighborhood  of  p^  and  a  nega¬ 
tive  slope  in  the  neighborhood  of  p_  .  While  shocks  w^th  initial 

Co 

pressure  higher  than  the  pressure  at  this  sonic  point  bul\ld  up  to  HVD, 
shocks  with  lower  pressures  decay  to  LVD.  In  most  cases,,  however,  HVD 
propagating  in  the  cavitated  liquid  will  eventually  becomes  HVD  propa¬ 
gating  in  virgin  liquid.  HVD  in  cavitated  liquid  will  usually  produce 
HVD  in  virgin  liquid  because  it  outruns  the  cavitation  field  and  produces 
a  shock  in  virgin  liquid  that  is  strong  enough  to  produce  HVD. 
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Construction  of  Hugoniot  Curves  for  the  Initiation  and  Propagation  of  LVD 


Relationships  for  Hugoniot  curves  with  a  positively  sloped  lower 
branch  and  a  negatively  sloped  upper  branch  were  fomulated  to  account 
for  the  initiation  and  propagation  of  LVD.  The  condition  for  a  Hugoniot 
to  have  a  positive  slope  will  be  derived  first,  however,  because  Hugoniot 
curves  are  usually  found  to  have  a  negative  slope.  The  equation  for  the 
slope  of  a  Hugoniot  is  readily  obtained  as 


dp  r/deN  .  dXn 

^  Lw) ,  -  *<vo  - v)  - 

v,  A 


(16) 


by  eliminating  de  from  the  differential  forms  obtained  by  differen¬ 
tiating  Eq.  (2)  and  the  e  =  e(p,v,X)  equation  of  state.  Since  we  are 

considering  a  case  in  which  (3e/3p)  .  >  0,  (3e/dv)  .  >  0,  and 

v,  X  p,  X 

2(3e/dp)  >  (v  -  v),  it  follows  from  Eq.  (16)  that  the  slope  of  a 

v,  A  o 

reactive  Hugoniot  is  positive  only  when  the  reaction  term  Q  dX/dp  is 
large  enough  to  make  the  term  in  brackets  on  the  right-hand  side  of 
Eqt  (16)  negative. 


The  equation  for  a  reactive  Hugoniot  curve  for  our  polytropic 
material  is  readily  obtained  as 


r  (K  +  1) 

PLV  7T7U  -  voJ 


2QX 


P„TV  ^ 

CJ  o 
(K  -  1) 


(17) 


by  eliminating  e  from  Eqs.  (2)  and  (7),  and  making  use  of  Eq.  (11). 

Setting  X  =  0  in  Eq.  (17)  gives  v  =  (K  -  1)/(K  +  l)v  for  a  strong 

o 

nonreactive  shock,  and  setting  X  =  1  in  Eq.  (17)  gives  the  equation 

for  the  fully  reacted  Hugoniot  curve.  The  condition  p  X*/p*  =  1 

CJ 

for  a  sonic  point  on  a  Hugoniot  curve  is  readily  obtained  by  setting 
v  =  tKvQ/(K  +  1)]  in  Eq.  (17).  It  is  convenient  for  computational 
purposes  to  set  K  =  3,  mid  rewrite  Eq..  (17)  as 
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X  ^CJ 

4  p 


(18) 


-  i  = 


Equation  (18)  gives  the  equation  for  the  slope  of  this  Hugoniot  as 


iE  (v  .  I»  . 

dv  \  2  dp/ 


US) 


and  it  follows  from  Eq.  (19)  that  the  Hugoniot  has  a  positive  or  a  nega¬ 
tive  slope  depending  on  whether  v  <  (v  /2  +  QdX/dp)  or  v  >  (v  /2  +  QdX/dp). 

o  o 

It  follows  from  Eq.  (18)  that  the  specification  of  either  a  v  =  v(A) 

relationship  or  a  p  =  p(A)  relationship  defines  a  reactive  Hugoniot 

curve.  In  the  present  work,  however,  X  =  A(p)  relationships  were 

formulated  to  generate  Hugoniots  with  the  properties  required  to  account 

for  the  initiation  and  propagation  of  LVD  and  for  its  transition  into 

HVD.  For  notational  simplicity,  v  and  p  _  were  used  to  denote  the 

o  CJ 

initial  specific  volume  of  the  cavitated  liquid  and  its  CJ  pressure, 

and  v  '  and  p '  were  used  to  denote  these  quantities  for  virgin 
o  CJ 

liquid.  Values  for  the  initial  volumes  of  virgin  and  cavitated  liquid, 

the  HVD  propagation  velocity,  and  the  critical  pressure  for  the  onset  of 

reaction  were  chosen-,  respectively,  as  v  r  =  1  cc/gm,  vq  =  1.25  cc/gm, 

U  =6  mm/fJsec,  and  p*  =  2  kbar.  The  CJ  condition  gives  c  =  4.5  mm/fjsec 
D  i  CJ 

and  u  =1.5  mm/JJsec,  and  Eq.  (11)  gives  p  =  90  kbar  and  p  = 

CJ  CJ  CJ 

72  kbar  for  the  CJ  pressure  attained  in  virgin  and  cavitated  liquid, 
respectively. 

The  unreacted  and  fully  reacted  Hugoniot  curves  for  cavitated  liquid 

obtained  by  setting  X  =  0  and  X  =  1  in  Eq.  (18)  are  shown  in  Figure  1. 

The  unreacted  curve  v/v  =  0.5  is  shown  as  NR,  and  the  fully  reacted 

o 

curve  as  ECD.  The  point  N  represents  the  criticality  condition  for 
onset  of  reaction  in  the  cavitation  field,  and  the:  point  C  represents 
the  CJ  point  for  HVD.  The  sonic  locus  v/vq  =  0.75  is  also  shown  in 
Figure  1  as  'SC. 
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PRESSURE  —  kbar 


i 


ytfibUi  //  ;;  /*  'hmtpff' 


n  r^I^Ltion^hlp  -f or  qoiv«.t3^tiqg.  r£%ef  iv#^ugc>iviofc  curve? 

f iS?s^  f Gnsuistted.  Bljsee  A  =  0  when,  p  -  -p*,- 
&gd:  X_=  -l  WMb_  J>  "  S(j3  p?r  these  ciirvcs,,  thrj  X  =  ACp)  relationship 
was^cHose^siffiply  afe  x_.  .  X 

S'  .  ..  F1  -  ^  T  j* 

*’■"  j  1  - !>*/>CJj  -;.'  _  <80> 

to-satl'gfj'^ihese  eojndiEions,  Vva;i?3es_qiv'thq^para3ietsr  n  are  restricted, 

liQwpyejrj^by  the  cajxdi.vJoss  iffiat/the ^slop'e^f  ^  reactive  Hugbaibi  carve 

-must  satisfy  at  -;e'  .a^!f .  K'.:  ^V:£3.'ues  pt  ii  .  must  he  chosen  so  that  the 

reactive  KugojiiDjt _  cu>.ye<s'afisii'es?the  ip^  condit  ion  (dp/dv)  =  -  S(p/v')- 

at  C.  ~  Thov:!nq#-ai<o^be':ehqsen.  so  ^li/DX  >  0  at  S  to  ensure  that 

shocks  satisfy  the/eohdiiion  Op/jDy  ^  0.  Equation  (19)  shows  that  a 

relative  Tiugoiriot:  satisfies  the  -C3  cchciltLid.’i  if-  dX/tip  =  p  when  X  =?  I;, 

and-  differentiating,  Sa^  020)  shows  that  the  letter,  condition  is  satisfied 

when  h  >  1  .  Di^-ferentiafcion  of  the.  equation  iluVv  =  <-p  r  Xp(,y'2)  obtained 

by  combining  Sqs.  (2)  and  (18)  gives  h  <2(1  -  p*/p„ )  ~  1.944  as  the 

condition  that  n  insist  -satisfy  to  ensure  that  Pu/Da  >0  at  H,. 

Values  of  n  k  2(1  -  p*/p  )  give  Dp/Du  <  0  at  II,  and  the  value 

i  CJV  -  - 

n  =  2,(1  r'.  pf/p  );  gives  a  Hugoniot  .in- the  (p-u)  plane  with  an  Infinite 
slope  at  p  =  p*V  The  additional  conditions  on  the  derivatives  at  C 
and  N  therefore  restrict  the  values  of  p  to  the  range  1  <  n  <  1.944. 


Equations  (18)  and  (20)  were  used  to  calculate  reactive  (p~v) 
Hugoniots  for.  different  values  of  n;  the  results  are  given  in  Table  I. 
The  limiting  values  n  =  1  and  n  =  1,944  were  used  to  determine  the 
domain  of  the  (p-v)  plane  covered  by  the  family  of  reactive  Hugcniots 
defined  by  Eq.  (20).  The  Hugoniot  curves  calculated  with  a  »  1,  1.3, 
and  1.944  are  shown  in  Figure  1.  Examination  of  Figure  1  shows  that 
while  the  curve  with  n  =  1  does  not  have  a  sonic  point  below  C,  the 
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Table  I 

REACTIVE  HUGONIOT  CURVES  FOR  CAVITATH)  LIQUID 


n  = 

1 

n  = 

1.3 

n  = 

1.5 

n  = 

1.944 

p 

(kbar) 

V 

(cc/gm  5 

X 

V 

(cc/gm) 

X 

V 

(cc/gm  ) 

X 

V 

(cc/gm) 

X 

3 

) 

0.732 

0.014 

0.764 

0.019 

0.786 

0.021 

0.832 

0.028 

4 

0.785 

0.029 

0.833 

0.037 

0.864 

0.042 

0.932 

0.055 

6 

0.839 

0.057 

0.901 

0.074 

0.942 

0.084 

1.030 

0,108 

8 

0.866 

0.086 

0.934 

0.110 

0.979 

0.126 

1.075 

0.160 

10 

0.882 

0.114 

0.954 

0.146 

0.100 

0.167 

1.098 

0.210 

30 

0.925 

0.400 

0.989 

0.485 

1.026 

0.535 

1,097 

0.630 

50 

0.9336 

0.686 

0.975 

0.778 

0.996 

0.824 

1.028 

0.895 

70 

0.937 

0.971 

0.943 

0.990 

0.345 

_ a, 

0.995 

0.946 

0.999 

curve  with  n  =  1.3  has  a  sonic  point  I  at  a  pressure  of  «  8  kb&r, 

/ 

and  the  curve  with  a  =  1*944  has  a  sonic  point.  I  at  a  pressure  of 
as  4.5  kbar.  Since  LVD  pressures  *Jk4%  Cue  6  to  9  kbar  range,  values  of 
n  in  the  range  1  <  n  <  1.944  can  be  used  to  model  LVD.  The  sonic 
point  1  (p  =  8.3  kbar,  X  =  0.115)  on  the  reactive  Hugoniot  with  n  =  1.3 
was  chosen  to  represent  LVD  because  this  state  is  attained  in  a  shock 
propagating  at  a  velocity  of  2.1  mm/fjsec  and  the  typical  propagation 
velocity  of  LVD  is  U*  «  2  rtm/psec. 

We  are  now  in  a  position  to  construct  Hugoniot  curves  with  either 
one  sonic  point  (pi ,  Xi)  or  with  two  sonic  points  (pi ,  Xi)  and  (p2,  X2) 
which  were  proposed  earlier  in  this  report  to  account  for  LVD.  The  case 
when  a  Hugoniot  has  one  sonic  point  will  be  considered  first. 

In  the  first  case,  pi  =  8.3  kbar  and  Xi  =  0.115  since  I  was 

chosen  to  represent  LVD.  The  lower,  reactive  branch  of  the  Hugoniot  with  a 

positive  slope  is  determined  by  Eqs.  (18)  and  (20),  with  n  =  1.3  and  with 

p  and  A  subject  to  the  conditions  (pi  =  2  kbar)  £  p  S  (p*  =8.3  kbar) 

and  0  £  X  £  0.115.  The  upper,  nonreactive  branch  of  the  Hugoniot  with 

* 

a  negative  slope  is  determined  by  Eq.  (18)  with  p  £  p^  =  8.3  kbar  and 
X  =  Xi  =  0,115.  The  Hugoniot  curve  calculated  according  to  this  pre¬ 
scription  Is  shown  as  NIF  in  Figure  2  where  the  Rayleigh  lines 
representing  LVD-  and  HVD  are  shown  respectively  as  01  and  0C 

In  the  second  case,  pi  =  8.3  kbar  and  X2  =  0.115.  A  reactive 
Hugoniot  that  connects  N  and  I  and  has  two  points  of  intersection 
with  the  sonic  locus  SC  will  now  be  constructed.  It  is  clear  from  the 
reactive  Hugoniot  curves  shown  in  Figure  1  that  a  X  =  X(p)  relationship 
producing  Hugoniot  curves  with  the  desired  properties  can  be  obtained  by 
scaling  Eq.  (20X.  A  a  =  X(p)  relationship  for  reactive  Hugoniot  curves 
with  two  sonic  points  was  therefore  formulated  by  setting 
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FIGURE  2  HUGONIOT  CURVES*  TO  ACCOUNT  FOR  THE  INITIATION  AND 
PROPAGATION  OF  liVD 
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X  = 


(21) 


Application  of  an  argument  similar  to  the  one  used  previously  for  curves 

with  one  sonic  point,  however,  shows  that  n  in  this  case  is  restricted 

$  * 

to  the  range  1  <  n  <  2(1  -  p  /pa)  =  1.52.  The  lower  reactive  branch 

j- 

of  the  Hugoniot  curve  with  two  sonic  points  is  determined  for  these  val¬ 
ues  of  n  by  Eq.  (18)  and  (21),  with  p  and  X  subject  to  the  condi¬ 
tions  (p*  =  2  kbar)  ^  p  ^  (p^  =8.3  kbar)  and  0  ^  X  s  0.115.  The 
upper  frozen  branch  of  the  Hugoniot  is  determined,  as  before,  by  Eq.  (18) 
with  p  ^  p^  =  8.3  kbar  and  X  =  Xa  =  0.115.  Reactive  Hugoniot  curves 
connecting  N  and  I  were  calculated,  and  the  results  of  the  calcula¬ 
tions  performed  for  n  =  1.3  and  1.5  are  shown  in  Table  II.  The  reac- 
tiire  curve  calculated  with  n  =  1.5,  is  shown  in  Figure  2  as  NTI ,  and 
the  corresponding  Hugoniot  constructed  to  account  for  the  initiation  and 
propagation  of  LVD  is  shown  as  NTIF.  The  lower  sonic  point  is  located 
at  the  point  T  with  a  pressure  p*  ~  4.5  kbar. 

Since  p*  ~  4.5  kbar  and  ps>  =  8.3  kbar,  LVD  modelled  with  NTIF 
will  exhibit  the  following  properties.  Shocks  with  initial  pressures 
p^  <~4.5  kbar  w;ill  decay,  shocks  with  initial  pressures  in  the  range 
8.3  >  p^  >  ~  4.5  kbar  will  build  up  to  LVD,  and  shocks  with  p^  >  8.3 
but  less  than  the)  critical  pressure  for  onset  of  bulk  reaction  will 
decav  to  LVD. 

Although  the:  transition  of  LVD  into  HVD  will  not  be  considered 
further  in  this  report,  it  is  clear  that  X  =  X(p)  relationships, 
producing  curves  similar  to  the  reactive  Hugoniots  shown  in  Figure  1, 
can  easily  be  formulated  to  account  for  this  type  of  behavior. 
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Table  II 

REACTIVE  HUGONIOT  CURVES  FOR  LVD 


p 

(kbar) 

n  = 

L.3 

n 

?  1.5 

V 

(cc/gm) 

X 

V 

(cc/gm) 

X 

1 

2 

0.625 

0.0 

0.625 

0.0 

3 

0.799  1 

0.023 

0.823 

0.026 

4 

0.879 

0.045 

0.907 

0.050 

5 

0.920 

0.066 

0.947 

0.072 

6 

0.940 

0.084 

0.963 

0.090 

7 

0.948 

0.100 

0.961 

0.105 

8 

0.943 

0,113 

0.946 

0.114 

0.938 


0.115 


0.938 


0.115 


EXPERIMENTAL  STUDY  OF  LVD 


The  objective  of  the  experimental  phase  of  the  program  was  to  obtain 
a  more  accurate  and  quantitative  description  of  LVD.  Shock  wave  experi¬ 
ments  using  multiple  ytterbium  stress  gages  were  designed  and  performed 
to  achieve  this  objective.  A  new  method  of  gage  emplacement  was  developed 
to  allow  the  gages  to  move  with  the  flow  so  that  the  experiments  would 
provide  stress  histories  at  different  Lagrangian  positions  in  a  liquid 
undergoing  LVD. 

The  new  method  of  gage  emplacement  is  a  definite  improvement  over 
12 

the  previous  method,  which  was  developed  for  AFOSR  under  Contract 
F44620-69-C-0070  and  used  to  study  LVD  in  ethyl  nitrate  and  FEFO.  In 
these  previous  experiments,  the  ytterbium  gages  were  placed  in  the 
direction  of  the  flow.  The  gage  records  were  suitable  for  measuring 
shock  arrival  time  and  .for  determining  propagation  velocity,  but  later 
interpretation  was  difficult  because  it  was  not  known  whether  the  gages 
recorded  predominantly  Lagrangian  or  Eulerian  stress  histories.  Further¬ 
more,  placement  of  gages  in  the  direction  of  the  flow  caused  their 
response  time  to  be  large  (>'250  nsec).  Consequently  in  the  current 
program,  the  gages  were  placed  perpendicular  to  the  direction  of  the 
:f low  so  that  they  would  move  with  the  liquid  and  record  Lagrangian  stress 
histories  with  a  smaller  response  time. 

Six  experiments  were  performed  in  the  SRI  explosives  vault  facility 
using  a  mixture  of  tetranitromethane  (4  parts  by  weight)  and  nitro- 
methane  (5  parts  by  weight)  in  a  steel  tube  into  which  several  ytterbium 
stress  gages  had  been  emplaced  to  record  Lagrangian  stress  histories. 


A  schematic  diagram  of  the  experiments  is  shown  in  Figure  3.  A  cylindrical 

seamless  cold-rolled  steel  tube  was  used  to  contain  the  liquid.  The  tube 

was  14-in.  long,  had  an  inside  diameter  of  3/4  in. ,  a  wall  thickness  of 

1/8  in.,  arid  was  filled  with  liquid  to  within  an  inch  of  its  top.  The 

explosive  donor  system  consisted  of  a  bridgewire  detonator,  a  tetryl 

booster  pellet,  a  2-in. -diameter  by  2-in. -thick  tetryl  explosive  pad, 

* 

and  a  3-in. -diameter  Homalite  attenuator.  The  detonation  in  the  tetryl 
pad  was  initiated  by  the  exploding  bridgewire  and  the  tetryl  pellet,  and 
the  thickness  of  the  attenuator  was  chosen  to  produce  the  desired  shock 
pressure  entering  the  liquid.. 

Two  to  four  stress  gages  were  positioned  from  2  to  8  inches  above 
the  bottom  of  the  tube.  A  typical  four-terminal  ytterbium  foil  grid  gage 
element  is  shown  in  Figure  3.  This  element  was  attached  to  a  thin  layer 
of  Homalite-impregnated  fiberglass  and  covered  with  a  thin  layer  of 
Homalite  to  form  a  total  gage  package  with  a  thickness  of  «  6  mils.  The 
gage  package  was  cut  to  fit  snugly  into  the  steel  tube,  with  two  segments 
removed  as  shown  in  Figure  3  so  that  the  tube  could  be  filled  with  liquid, 
and  was  supported  by  a  thin  layer  of  glue.  A  1/8-in.  hole  was  drilled 
through  the  tube  approximately  an  inch  above  each  gage  package  location, 
and  the  leads  from  each  gage  were  brought  out  through  this  hole.  This 
method  of  emplacement  ensured  that  the  gage  would  survive  long  enough  to 
record  the  LVD  pulse.  Had  vrie  gage  leads  been  brought  straight  out  at 
the  same  level  as  the  gage  package,  they  would  have  sheared  off  as  the 
liquid  flowed  along  the  inside  wall  of  the  tube. 

An  additional  ytterbium  stress  gage  was  placed  at  or  near  the 
Homalite  attenuator-liquid  interface  to  record  the  stress  pulse  entering 
the  liquid. 

*Trade  name,  Homalite  Corporation. 
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The  gage  records  were  calibrated  using  the  formula  given  by 
13 

Ginsberg  for  calibration  of  ytterbium  foil  for  the  region  from  0  to 
10  kbar. 

Of  the  six  experiments  attempted,  all  but  one  yielded  satisfactory 
data.  The  one  that  failed  was  shot  2,  using  a  2-in. -thick  Homalite 
attenuator.  Post-shot  recovery  of  the  tube  fragments  indicated  that  shot 
2  experienced  a  high-velocity  detonation  (HVD),  Two  other  shots  in  this 
series  used  the  same  attenuator  thickness  and  therefore  had  nearly  the 
same  stress  pulse  magnitude  enter  the  liquid;  these  two  shots  exhibited 
LVD. 

The  relevant  parameters  fcr  the  five  successful  shots  are  shown  in 
Table  III,  along  with  the  peak  stresses  recorded  by  the  various  gages 
and  the  velocity  of  the  compressional  part  of  the  LVD  wave.  The  five 
experiments  span  a  range  of  peak  stress  incident  upon  the  liquid  from 
3.6  to  16  kbar. 

The  oscillographs  from  the  stress  gages  in  shot  3  are  shown  in 
Figure  4.  Three  gages  recorded  in  the  liquid  at  2,  4,  and  6  inches 
above  the  liquid-Homalite  interface,  and  a  fourth  gage  recorded  just 
below  the  interface.  The  fourth  gage,  Figure  4(d),  shows  a  peak  stress 
equal  to  the  stress  in  the  Homalite  and  then  a  decrease  to  the  stress 
entering  the  liquid  (the  liquid  has  a  lower  shock  impedance  than  the 
Homalite,  hence  the  difference). 

All  the  in-liquid  stress  histories  recorded  for  all  five  shots  show 
features  similar  to  those  shown  in  Figure  4  (a),  (b),  and  (c).  These 
include  a  rapid  rise  or  shock  to  a  peak  pressure  of  from  6.0  to  9.1  kbar, 
a  shock  velocity  between  successive  gages  of  from  1^87  to  1*99  im/^lsec, 
and  a  decay  to  nearly  zero  stress  of  several  microseconds  duration.  There 
is  no  noticeable  variation  (to  within  the  statistics  of  these  experiments) 
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Table  III 


EXPERIMENTAL  PARAMETERS  AND  RESULTS  FOR  LVD  SHOTS 


Shot  No. 

1  3 

4 

5 

6 

Attenuator 
thickness  (in.) 

2  ,  3--1/32 

3-1/32 

2-1/32 

5-1/32 

Peak  stress  entering 
liquid3  (kbar) 

16  '6;1 

8.r4 

b 

3.6 

First  in- liquid  gage; 
Height  in  tube  (-in.)’ 

2  2'  . 

7. 

2 

2 

Recorded*  peak  stress 
{kbar)  - 

7.25  7^4  ' 

_e.i5 

7.2 

7.25 

Second  in-liquid  gage; 
Height  in  tube(in.). 

5-  '  -  .  4 

■4- 

4 

4 

Recorded  peak  ^stress 
jCkbar)  '  -  , 

8.15  $.6 

7.5 

6.0 

10. 3C 

"  dv  -  - 

Measured  velocity  - 
■Om/'fJs'ec), 

lijo-  ;a*35  ' 

1,96 

1.99 

1.96 

Thlrcrih-1 igutd-  gage;  ;V- 
Height,  in:  tube  (in.  >  : 

:  ~  -  ,6 

6 

— 

6 

Recorded  peak  stress, 
(kbar) 

*•”*-  ;9'.  1 

8^25 

— 

8.6 

Measured  velocity  ' 
(mm/fteec) 

—  1.93 

1.87 

— 

1.92 

Fourth'  in-liquid  -.gagers 
Height  in  tube  (ih«)  — 

--  '  9 

8 

8 

Recorded  peak  stress 
(kbar) 

--  :  b 

7.8 

8.15 

Measured  velocity 
(mm/^sec) 

1.32  - 

;  1,91 

-- 

1.89 

Average  peak  stress  for 
iti-liquld  gages  (kbar) 

7:.r  7.7  1 

..  3  0 

6.6 

8.0 

Average  measured  velocity 
(mm/jUsec) 

11,96  1.93 

'  1.91 

1.99 

1.92 

a.  As  measured  by  gage  at  HomalitS"liqutd  interface. 

b.  Data  Of  insufficient  quality  to  axialyze^ 

c.  Questionable  data*  not  used  in  'averages « 

d«  Measured^  velocity  is  that  of  the  midpoint  of  the  compressive 
shock  between  this  and  the  previous  gage  location. 


M 

.  « 

l  ■ 
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in  peak  pressure  or  shock  velocity  as  a  function  of  incident  stress  into 
the  liquid.  In  other  words  the  LVD  wave  attains  the  sane  approximate 
peak  pressure  and  shock  velocity  regardless  of  the  magnitude  of  the 
initiating  shock.  Furthermore  this  approximately  steady  wave  is  nearly 
attained  within  the  first  two  inches  of  the  liquid,  as  shown  by  the 
closeness  of  the  average  peak  stress  measured  at  the  2-in.  location 
(7.4  kbar)  compared  with  the  overall  average  peak  stress  (7.6  kbar)  and 
the  similarity  of  the  average  shock  velocity  measured  between  the  2-in. 
and  the  next  gage  location  (1.96  mm/psec)  to  the  overall  average 
(-1.93  mm/psec). 

Figure  5  is  a  plot  of  the  arrival  times  of  the  midpoint  of  the  LVD 
shock,  in  microseconds  from  the  time  the  shock  enters  the  liquid,  as  a 
function  of  the  Lagrangian  position  or  initial  location  of  the  gage  plane 
in  the  liquid.  The  lines  for  all  five  shots  are  very  straight  and  nearly 
parallel,  indicating  constant  velocity.  The  only  difference  is  in  their 
vertical  position  on  the  graph.  For  example,  the  LVD  wave  for  shot  6 
arrives  3  psec  later  than  that  for  shot  1.  Since  the  stress  entering 
the  liquid  in  shot  1  is  16  kbar,  while  that  in  shot  6  is  only  3.6  kbar, 
the  difference  in  arrival  times  is  caused  by  a  combination  of  the  fol¬ 
lowing  two  factors:  (1)  The  shock  velocity  of  the  unreacted  liquid 
before  LVD  initiation  is  higher  at  higher  stresses  because  of  the  con¬ 
cave  upward  nature  of  the  liquid  Hugoniot.  (2)  The  LVD  might  take 
longer  to  initiate  at  lower  stresses.  Further  experimentation  is 
necessary,  however,  to  determine  which  of  these  factors  is  more 
important. 

The  in-\iquid  stress  gage  records  for  shots  3  and  6  were  digitized 
and  plotted  together  in  Figures  6(a)  and  6(b),  respectively.  The  hori¬ 
zontal  or  time  axis  was  shifted  so  that  th(  first  motion  or  the  foot  of 
the  compressional  shock  for  all  three  gages  coincided  at  zero  time. 
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DISTANCE  ALONG  TUBE  AXIS  —  inches 

MA-2383-3 


FIGURE  5  DISTANCE-TIME  PLOT  OF  MIDPOINT  OF  LVD  COMPRESSIONAL  WAVE 
FOR  ALL  5  EXPERIMENTS 
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(b)'FROM  SHOT  6 

MA-2383-3 

FIGURE  8  COMPOSITE  STRESS  HISTORIES  OF  IN-LIQUID  YTTERBIUM 
GAGES  (concluded) 
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If  LVD  were  precisely  steady  state  throughout  the  region  of  measurement, 
each  point  in  the  wave  would  travel  at  the  same  velocity,  and  as  a 
consequence  the  recorded  stress  histories  would  be  exactly  congruent. 

The  actual  results  are  fairly  congruent  in  each  case  except  in  two 
regions — near  the  peak  and  at  the  final  third  of  unloading.  The  peak 
stress  of  the  first  gage  is  significantly  lower  than  that  of  the  other 
two  for  both  shots  3  and  6,  indicating  that  perhaps  the  LVD  has  not  yet 
reached  its  maximum  or  steady-state  stress  level  in  the  first  two  inches 
of  liquid.  The  unloading  is  slower  for  gages  farther  up  the  tube,  which 
indicates  that  the  release  wave  velocity  is  less  than  that  of  the  com- 
pressional  wave.  The  lower  velocity  of  the  release  wave  is  qiii'te  likely 
because  the  liquid  is  in  a  heated  state  following  detonation,;  however, 
the  difference  is  very  small.  In  shot,  3,  the  release  velocity  at  4  kbar 
is  only  1%  lower  than  the  shock  velocity,  arid  from  there;  the  difference 
increases  to  a  maximum  of  8%  at  1  kbar.  For  shot  6,  the  differences  are, 
even  smaller. 

In  conclusion,  we  have  successfully  obtained  Lagrangiaii  stress 
histories  in  a  liquid  undergoing  LVD;.  The  records  indicate  that  the 
compressional  part  of  the  LVD,  wav.e  propagates  at  very;  nearly  constant 
peak  stress  and  velocity,  whfle  the  release  portion;  of  the  LVD  wave, 
although  initially  nearly  steady  state,  begins  to  /diverge  (i.e.,  propa¬ 
gate  at  slightly  lower  velocities);  as  the  stress  decays  to  zero.  This 
technique  can  now  be  used  with  a  high  degree  of  certainty  of  success  to 
study  in  detail  the  initiation  and  propagation  of  LVD. 


CONCLUSIONS 


Considerable  progress  has  been  made  in  a  basic  understanding  of  LVD 
and  in  determining  conditions  for  its  initiation  and  propagation. 

The  treatment  of  reactive  shocks  formulated  in  the  theoretical  study 

provides  for  the  first  time  a  simple  model  to  account  for  the  initiation 

and  propagation  of  LVD.  Shocked  states  are  described  by  Hugoniot  curves 

that  exhibit  two  sonic  points  in  the  region  of  2  to  10  kbar.  The  sonic 

♦  _ 

point  with  \the  higher  pressure  p2  represents  LVD,  and  the  sonic  point 
with  the  low;er  pressure  p*  defines  the  initial  pressure  needed  for 
shocks  with  a^positive  particle  velocity  gradient  to  build  up  to  LVD. 

If  pi  denotes,  the  initial  shock  pressure *  then  shocks  decay  when 
p^  <  P* i  build  v'ip  to  LVD  when  p*  s  p^  <  pa,  and  decay  to  LVD  when 
Pt  >  P2. 

The  development  of  multiple  ytterbium  gages  to  determine  Lagrangian 
stress  histories  in  shocked  liquids  provides  a  means  of  obtaining  a  more 
accurate  and  quantitative  description  of  the  states  attained  during  the 
initiation  and  propagation  of  LVD.  The  stress  histories  recorded  by 
the  ytterbium  gages  in  a  tetranitromethane-nitromethane  mixture  lead  to 
the  following  conclusions  about  LVD  in  this  liquid.  The  compressional 
part  of  the  LVD  wave  propagates  at  constant  velocity  with  a  constant 
peak  pressure;  the  flow  in  the  release  part  of  the  wave  is  steady  state 
as  the  pressure  falls  to  one-half  or  one-third  of  the  peak  pressure  and 
then  becomes  unsteady. 

The  results  of  the  experiments  show  that  LVD  in  the  tetranitromethane- 
nitromethane  mixture  is  initiated  by  shocks  with  initial  pressures  the 
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3.6  to  16  kba?  range,  and  the  reactive  Hugoniot  curves  constructed  in 
the  the  theoretical  study  lead  to  the  conclusion  that  shocks  will  build 
up  to  LVD  when  their  initial  pressures  are  greater  than  ~  4.5  kbar,  and 
will  decay  to  LVD  when  their  initial  pressures  are  greater  than  ~  8.5  kbar. 

Additional  work  is  required,  however,  to  develop  a  more  fundamental 
treatment  of  LVD  because  the  work  on  reactive  shock  discontinuities  is 
semiempirical ,  and  because  the  computer  calculations  on  burning  bubbles 
did  not  provide  a  condition  for  calculating  the  degree  of  incomplete 
reaction  in  the  compressed  cavitation  field.  To  be  more  specific,  LVD 
pressure  is  a  parameter  that  must  be  known,  and  the  critical  pressure 
for  onset  of  reaction  in  the  cavitation  field  and  the  relationship  for 
the  reactive  Hugoniot  curve  were  chosen  rather  than  calculated  from  first 
principles.  The.  calculations  on  burning  bubbles  did  not  provide  condi¬ 
tions  for  evaluating  these  quantities  because  the  calculations  were 
performed  for  a  contracting  bubble  in  an  incompressible  liquid,  and  the 
reaction  is  probably  extinguished  during  bubble  expansion  when  compressi¬ 
bility  factors  cannot  be  ignored.  It  is  therefore  necessary  to  formulate 
a  satisfactory  model  for  determining  the  amount  of  reaction  in  the  com¬ 
pressed-  cavitation  field  as  a  function  of  pressure  so  that  all  these 
quantities  can  be  calculated.  A  study  of  reactive  bubbles  in  a  compress¬ 
ible  liquid  undergoing  contraction  and  expansion  should  be  continued  to 
develop  such  a  model.  A  study  of  reactive  shock  propagation  including 
divergence  and  reaction  behind  the  shock  should  also  be  undertaken  to 
test  the  assumption  that  LVD  can  be  treated  as  a  one-dimensional  wave. 
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Appendix 


CALCULATION  OF  THE  COLLAPSE  OF  A  VAPOR  BUBBLE 
LEADING  TO  HOT-SPOT  INITIATION 

Introduction 

The  important  role  of  bubble  hot-spot  ignition  in  the  initiation  and  pro¬ 
pagation  of  low  velocity  detonation  (LVD)  was  discussed  in  a  previous  report,  and 
equations  were  derived  that  govern  the  collapse  of  a  spherical  vapor  bubble.*  Of 
interest  here  is  the  collapse  of  a  bubble  leading  to  high  interior  temperatures  caus¬ 
ing  local  initiation  .of  reaction  which  may  ultimately  Ignite  the  bulk  material  lead¬ 
ing  to  explosion. 

Theoretical  studies  of  hot-spot  initiation  that  are  related  to  the  adiabatic 
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compression  of  a  bubble  have  been  carried  out  by  Zinn ,  Enlg,  and  Gill.  These 
'  are  computations  of  ignition  delay  time  which  are  applicable  to  a  bubble  with  a  sta¬ 
tionary  radius  and  with  no  mass  transfer  at  the  bubble  wall.  In  the  present  compu¬ 
tations,  however,  we  are  Interested  in  the  temperature  rise  and  the  onset  of  reaction 
or  ignition  in  a  collapsing  vapcr  bubble.  We  include  here  the  effects  of  both  heat 
and  mass  transfer  (vaporization)  at  the  bubble  wall,  and  formulate  the  problem  in 
Lagrangian  coordinates  that  move  with  the  bubble  wall  during  collapse.  It  is  assumed 
that  both  the  vapor  and  liquid  are  nonvlscous  and,  prior  to  collapse,  they  have  the 
same  temperature  Tq.  The  interior  of  the  vapor  bubble  is  assumed  to  be  uniform, 
an  assumption  that  is  justified  provided  that  the  thermal  diffusivlty  of  the  vapor  is 
much  greater  than  that  of  the  liquid  and  that  the  speed  of  the  bubble  wall  motion 
is  much  less  than  the  sound  speed  in  the  Vapor.  The  first  condition  is  certainly 
true  for  most  substances  for,  although  the  thermal  conductivity  of  the  liquid  is  a- 
bout  ten  times  greater  than  that  of  the  vapor,  the  thermal  diffusivlty  of  the  vapor  is 
still  much  greater  than  that  of  the  liquid.  The  second  condition  has  been  discussed 
by  Hickling5  in  connection  with  the  collapse  of  a  bubble  of  inert  gas  in  a  liquid. 
Hickling  found  that  the  sound  speed  in  the  gas  was  indeed  much  higher  than  the 
velocity  of  the  bubble  wall  for  initial  bubble  radii  of  Rq  =  10" 1  to  10~3  cm.  and 
overpressures  of  about  4  atm.  Under  such  conditions,  the  pressure  in  the  vapor 
bubble  remains  approximately  uniform  because  the  effects  of  disturbances  at  the 
wall  will  have  time  to  be  propagated  throughout  the  bubble  interior.  It  is  clear 
that  during  the  early  stages  of  collapse  the  assumption  of  uniform  bubble  Interior 
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will  certainly  be  valid  in  the  present  calculation.  Further,  to  keep  the  problem 
tractable,  a  spherical  bubble  will  be  assumed  even  though  shock  compression  in 
general  produces  asymmetric  collapse  and,  in  addition,  it  s  known  that  spherical 
bubbles  are  generally  unstable  during  the  latter  stages  of  collapse. 

The  formulation  of  the  equations  to  be  solved  is  given  in  the  following 
section,  followed  by  a  discussion  of  the  numerical  methods  used  to  solve  the 
equations.  Th?  results  of  several  example  solutions  are  then  discussed. 
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Formulation  of  the  Problem 


A  very  general  formulation  of  the  equations  governing  the  motion  of  a 
bubble  in  a  liquid  was  given  by  Hsieh,  as  well  as  simplified  equations  for  a 
spherical  bubble  with  a  uniform  interior.  With  the  above  assumptions  and  the  ad¬ 
ditional  assumptions  that: 

(1)  No  body  force  is  present, 

(2)  There  is  no  translational  motion  of  the  bubble  relative  to  the  liquid;, 

(3)  Viscosity  coefficients  are  zero  throughout; 

(4)  The  liquid  is  incompressible; 

(5)  The  equation  for  the  temperature  in  the  liquid  T^  (r,  t)  can  be  written 


9Ti  .  „  ,  2  9Te 

T5T  +  vi  "ST  =  “i  \"S?"  +  T  "ST, 


(1) 


where  r  is  the  radial  coordinate  from  the  bubble  center  and  a ^  is  the  thermal  dif- 
fusivity  of  the  liquid,  which  will  be  assumed  constant.  Since  the  liquid  is  treated 
as  incompressible,  and  the  density  of  the  vapor  pv  is  generally  much  less  than  that 
of  the/liquid,  from  the  continuity  requirement  the  liquid  velocity  may  be  directly 
related  to  the  bubble  radius  R(t)  as  follows: 


(2) 


lhat  is,  the  contribution  of  mass  transfer  (vaporization)  at  the  nubble  surface  to  R 
is  neglected.  Note  that  the  assumption  of  uniform  bubble  interior  requires  that  the 
quantities  inside  the  bubble  pv,  Tv,  Pv  are  functions  of  time  only.  The  initial  con¬ 
dition  will  be  that  at  t=0  a  bubble  of  radius  RQ  Is  at  rest,  with  the  vapor  inside  at 
initial  density  and  pressure  of  pvQ,  respectively.  The  assumption  of  uniform 
initial  temperature  Tq  for  the  system  requires  that 

Tv(°)  =  tyr,0)  =  Tq  .  (3) 


Further,  the  boundary  condition  at  the  bubble  wall  requires  that 


Ty(t)  =  T^(0,t) 
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as 


and 


when  L  is  the  latent  heat  per  unit  mass#  is  the  specific  heat  at  constant  volume 
of  the  vapor,  and  X^  is  the  thermal  conductivity  of  the  liquid.  The  quantities  X^, 

L,  and  Cv  are  assumed  constant,  as  are  Q,  Z,  and  E  which  denote,  respectively, 
the  specific  heat  of  reaction,  the  frequency  factor,  and  the  activation  energy.  The 
last  term  in  Eq.  (5)  is  the  heat  release  in  the  bubble  due  to  gas  phase  reaction  which 
is  generally  negligible  until  Tv  reaches  the  order  of  E/fe. 

Since  the  (uniform)  temperature  in  the  vapor  bubble  is  determined  by  the 
temperature  of  the  liquid  at  the  bubble  wall,  it  follows  that  the  (uniform)  pressure 
in  the  bubble  ?v  is  given  by  the  vapor  pressure  of  the  liquid  at  the  wall  or  bubble 
temperature.  We  will  therefore  make  the  additional  assumption  that  the  vapor  pres¬ 
sure  is  the  equilibrium  vapor  pressure  at  the  bubble  wall  and  that  the  bubble  con¬ 
tains  pure  vapor.  Thus,  the  possibility  that  absorbed  gas  in  the  liquid  may  con¬ 
tribute  to  the  interior  pressure  of  the  bubble  is  excluded,  and  the  results  will  be 
applicable  to  a  pure  liquid  or  to  cases  in  which  the  partial  pressure  of  absorbed 
gases  is  negligible  in  comparison  with  the  vapor  pressure.  The  equilibrium  vapor 
pressure  requirement  is  certainly  valid  in  the  early  stages  of  collapse  until  the 
rate  of  change  of  the  pressure  in  the  bubble  becomes  comparable  to  the  kinetic 
rate  of  vaporization.  To  enable  us  to  treat  the  problem  analytically,  we  will  relate 
the  vapor  pressure  of  the  liquid,  .md  hence  the  pressure  in  the  bubble,  to  the 
temperature  at  the  bubble  surface  by  use  of  the  Calsuis-Clapeyron  equation 


P 


v 


-P0exp 


(6) 


where  ft  is  the  ideal  gas  constant  and  PQ  is  the  vapor  pressure  at  the  initial  temp¬ 
erature  T  .  Further,  the  thermodynamic  quantities  in  the  vapor  phase  will  be  re- 

o  % 

lated  by  the  ideal  gas  equation  of  state 

P„  -  *  Pv  Tv  .  (7) 
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The  additional  equation  required  to  complete  the  formulation  of  the 
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problem  is  the  differential  equation  for  the  motion  of  the  bubble  wall 


3  .  „  (P.-pJ 

R  R  +  4  (R)  = - -  ,  t >0, 

Z  Pjt 


(8) 


where  ?m  ,  assumed  constant,  is  the  Initial  pressure  in  the  bulk  liquid  and  corres¬ 
ponds  to  the  pressure  at  large  distance  from  the  bubble  during  collapse.  The  in¬ 
itial  conditions  for  Eq.  (8)  are  taken  as 


and 


R(0)  =  R, 


R(0)  =  0 


(9) 


Equations  (1),  (2),  (6),  (7),  and  (8)  now  form  a  set  for  the  dependent 
variables  V^(r,t),  T^(r,t),  Pv(t),  and  py(t),  subject  to  the  conditions  of  (3),  (4), 
(5) ,  and  (9), 


Introduce  the  dimensionless  variables 


x  = 


gipr  (r3  -  R3  (t)) 


/V 

R 


=  ^c 
»  R(t)/R 


do) 


e„  = 


T  -  T 
v  ,  o 


9  = 


T,  -  T 
l  o 


and  the  parameters 


rj  =  L/RTq 

To  -  Vo /PJ1/2 

Th = 

rr  *  »T0/QZ, 
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where  r  is  the  characteristic  time  for  bubble  collapses  is  the  characteristic 
c  n 

time  for  heat  transfer  to  the  bubble  wall ,  and  rr  is  the  dimensionless  frequency 
factor  for  the  reaction  term  in  Eq.  (5). 

Equation  (1)  can  now  be  written  as 


H) 


and  the  conditions  (3),  (4)  ,  and  (5)  become  respectively. 


9V(0)  =  0(x,O)  =  0,x^0 


8„(t)  «  e(o,t)  ,  tao 


88  fn  ti  -  n  Cy  /pv°\  /Th  \  /  1  \  n  t  \ 

-jj  (0,t)  -  (v- 1)  j  \-g-j exp  7]  ^T+rJJ  x 

IdR  /_2_  \  1  /_1_\  ^v_  \(_JL-'\  ._L  +  _2_ 

,df  \1+ev  7  3  V1  +  9v/  dt  lA1+evj  ^  ’"1 

*  W(£)(“*  ['  A  i^k)}  6XP  KJ) 


where  y  is  the  ratio  at  specific  heats  of  the  vapor, phase  and  ovQ  is  the  initial  den¬ 
sity  in  the  vapor  phase. 


The  bubble  equation  (8)  can  be  written  as 


with  the  initial  conditions 


(#)  *  4  (■ 


R(0)  =  1 


\3 
dR  \ 


P  -  P 


(0)  =  0  . 
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uatlons  To  Be  Numerically  Integrated 

The  nondimensional  heat  conduction  equation  (12)  in  the  liquid  may  be 
written  as 


Two  boundary  conditions,  an  initial  condition,  and  an  equation  for  the  bubble 
radius  must  be  given  to  complete  the  specification  of  the  problem.  Far  away  from 
the  bubble  surface,  the  liquid  temperature  is  assumed  to  be  constant,  hence, 

0(x  -*  -)  =  0  (20) 

At  the  bubble  surface  the  energy  balance  (15)  can  be  written, 

aS. |  -  Vo  fJ!L)1//2  1  II  If  \  .  ,)  dR_ 

Sx  lx=0  “  XTq  V  5  V,™! 


+  T(I+^}  ^  '(l4%]-}  +  FT 

-  t  (_  mt"'  \exp  [-  $r~  (T+e7} 


¥~  -  exp 


where  the  pressure  ratio  Pj/Pq  Is  given  by 


T^=exp  =  (22i 

and  where  the  subscript  1  refers  to  conditions  at  the  bubble  wall.  The  initial 
condition  for  equation  (19)  is 


«£=  0)  =  0 


The  equation  for  the  motion  of  the  bubble  wall  is 

r4  +i<4){,+-£  4* -(!-.£-> 


2  at 


A  A 

R  dt 
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where, 

R(i=0)  =  l  %No)«a 

dt 


(25) 


Transformed  Equations 

The  temperature  gradient  is  greatest  at  the  bubble  surface  and  decays 
rapidly  with  distance  into  the  liquid.  For  this  reason  it  is  desirable  to  transform 
to  a  coordinate  system  in  which  constant  increments  in  the  transformed  plane 
correspond  to  very  small  increments  near  the  surface  and  ever  larger  ones  as  the 
distance  from  the  surface  increases.  Also,  the  boundary  condition  (20)  need  not 
be  applied  at  x  =  •,  since  during  the  very  rapid  bubble  collapse,  the  heat  trans¬ 
ferred  from  the.  bubble  interior  can  only  propagate  a  short  distance  into  the  liquid. 
There  are  possible  coordinate  transformations  possessing  the  desired  char¬ 
acteristics,  the  one-selected  for  use  in  the  current  analysis  is  given  by, 

1  +  if;  (ea  - » 


or. 


x 


qZ/aZ 


This  transformation  maps  the  region 


(26) 


x  =  0  •+ 


xmax 


into 


Z  =  0  -»  1 


where. 


xmax 


-p-  - 1) 

(ea  - 1) 


(28) 
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The  step  size  in  the  transformed  plane,  AZ,  is  constant.  In  the  physical  plane 

the  first  mesh  point  is  a  distance  Ax  from  the  bubble  surface.  The  ax  intervals 

o 

then  grow  in  an  exponential  manner  as  one  proceeds  away  from  the  surface.  The 
derivatives  of  the  transformation  are  also  required. 


(ea  -  X) 


then. 


A7 

uZ  rji  _  6 

■a?  *z  --^r 

j3  n  ~2j3Z 

0=z’"-w---z^ 


The  following  relations  can  then  be  used  to  transform  the  heat  conduction  equa¬ 
tion  (19). 

a  _7,  a 

(32) 

*  _  7„  a  .  7,»  a* 

9?-z  s'?  +z  aZ*- 


The  transformed  equation  becomes, 


M  i  (£  +  3X)1/3  [(£  +  3x)  Z"  +  4Z’j  ||-  +  (&  +  3x)  1 

St  n 


The  transformed  boundary  conditions  are: 


0  =  0  at  Z  =  1 


00 1  -  1  ®0  I 

"52'Z»0  Z’(6)  ^x"lx=0 


The  bubble  equations  (24)  and  (25)  are  not  affected  by  the  transformation. 
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Finite  Difference  Solution 

Equation  (33)  subject  to  the  boundary  conditions  (34)  and  (35).  and  the 
initial  condition  (23)  ,1?  solved  using  the  Crank-Nlcolson  Implicit  finite  difference 
scheme.  A  network  of  mesh  points  is  created  as  shown  in  Figure  1. 


►  t 


#  Known  Point 
O  Unknown  Point 

X  Equation  Evaluated 


Figure  1.  Finite  Difference  Mesh 


With  the  Crank-Nlcholson  method  the  derivatives  and  coefficients  in  the  differential 
equation  are  evaluated  at  the  midpoint  m+l/2,  n+1/2.  The  resulting  solutions 
are  accurate  to  second  order.  For  nonlinear  problems  the  difference  equations 
are  normally  linearized,  enabling  a  solution  to  be  achieved  by  solving  a  coupleu 
set  of  linear,  rather  than  nonlinear,  algebraic  equations.  In  such  cases,  one 
iteration  must  be  performed  to  regain  second  order  accuracy.  The  following  repre¬ 
sentations  of  the  derivatives  are  used: 


!/2 


3*6  _  I  /« 

w~ 1/2 


^m,n+l~  0m,h-l^  .  ^m+l,n+l  jn+I,n-l^ 
2A&  '  2AZ  ; 

(em  -  20  +  0m  .) 

>_r 


(0m+l.n+r  20m+l.n+  0m+l.n-l) 

— ^ - Hr1 - — 


(36) 


(37) 


50  0m+l,n  0m,n 

at  At 


(38) 
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After  replacing  all  of  the  derivatives  in  equation  (33)  by  their  central  difference 
analogs,  (36)  to  (38),  the  following  set  of  difference  equations  is  obtained. 


where 


An  ®m+l,n-l  +  Bn  0m+l,n  +  Cn  9m+l,n+l  =  Dn 


(39) 


n  =  2,3,  ...N-l 


An  =  (bn  '  en> 
Cn=  ~  (bn  +  en) 


1  =  surface 
N=  far  boundary 


B  *  i  +  2e 
n  n 


(40) 

D_  =  8  +b  0-  +ejr 
n  m,n  n  *  n  ZE  _ 
m,n  m,n 


%,n=9m'n+1"3m'n-1 


0V-7  -8  „ . ,  -  20  +0  . 

ZZm  n  m»n+l  m,n  m,n-l 


(41) 


and. 


T  (R®  +  3x) 
b  =  A  t  _  Av - n 


1/3  r 


n 


n 


+3xJ 


1/3 


e_  =  A  t 
n  r 


c  av  n 


n 


7Z7T 


(R*  v  +  3*J  K  +  4Z* 

av  n  n  n 


(K’av  +  3 V  Z’a 


(42) 


Once  the  boundary  conditions  have  been  applied,  (39)  yields  a  set  of  N-2,  tridia¬ 
gonal,  linear  algebraic  equations,  which  may  be  efficiently  solved  using  a  simple 
Gaussian  elimination  algorithm.  To  reduce  equations  (39)  to  tridiagonal  form  the 
boundary  conditions  at  Z=0  and  Z«1  must  be  used  to  evaluate  the  terms 

Vm+l,labdCN-iem+l,N  (43> 

The  far  boundary  condition  (34)  given  9m+1  n”0,  sur*ace  boundary  condition 
is  nonlinear  and  may  be  treated  in  numerous  ways.  The  method  currently  being 
used  Iteratively  selects  values  for  0m+j  ^  unt^  the  derivative  at  the  surface  given 
by  (35)  is  equal  to  the  calculated  derivative,  to  within  a  specified  tolerance.  The 

calculated  surface  derivative  is  evaluated  using  a  four  point  formula 

A-13 


(44) 


Mi 

8Z‘Z=0 


(-119rn+l.l  +  188m+l. 


2  ~  "m-Hf3  +  29m+l,4^ 


6AZ 


The  method  of  solution,  as  given,  assumes  that  the  instantaneous  bubble  radius 
is  known,.  In  reality,  however,  the  bubble  radius  is  a  function  of  pressure,  and, 
hence,  by(22),a  function  of  surface  temperature.  In  order -to  avoid  a  simultaneous, 
iterative,  solution  of  the  heat  conduction  equation  and  the  bubble  equation  (24), the 
calculation  of  the  bubble  radius  is  done  in  an  uncoupled  manner.  This  procedure 
for  uncoupling  the  bubble  radius  calculation  appears  to  work  satisfactorily,  and 
is  described  below. 

The  second  order  bubble  differential  equation  (24)is  rewritten  as  a  set  of 
two  first  order  equations, 

rM  +  -|.r2  +  .M.r=- (i  _  £L)  (45) 

dt  1  R  *• 

~  =  R  (46) 

dt 


a  simple  Euler  integration  of  (45)  and  (46)  is  then  carried  out. 


Rt+4'f  " 


A 

At 


(47) 

(48) 


When  calculated  in  this  manner,  the  bubble  radius  and  velocity  at  t+A  t  depend  only 
on  previously  known  values  at  time  =<t.  The  bubble  parameters  given  by  (47)  and 
(48)  are  then  used  in  evaluating  the  term  Rgv  in  (42),  and  the  surface  heat  transfer 
rate  (21). 
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Rq  =  ,1,  .OX  cm 

T  =  300°K 

o 

=  1.13  gm/cm3 

Pw  =  107,  10s  dynes/cm2 

ct  =  .45  cal/gm°K  =  1.9  x  107  erg/gm°K 

-  135  cal/gm  =  5.65  x  109  erg/gm 

M  =61  g/g-mole 
H  =  .0062  gm/cm  sec  (poise) 

X  =  2.15  x  104  erg/cm  sec°K 

E  =  2.25  x  1012  erg/mole,  or  1.125  x  1012  erg/mole 

=  3.69  x  10  11  erg/gm,  or  1.845  x  10u  erg/gm 

Z  =  3.98  x  1014  sec  “1 

Q  =  4.6  x  lO10  erg/gm,  or  0.0 

y  *  1.2 

B  =  8.317  x  107  ergs/g-mole0C 
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Example  Calculations  and  Results 


The  equations  tor  the  collapsing  spherical  vapor  bubble  were  numerically 
integrated  for  two  values  of  initial  bubble  radius  and  two  overpressures,  using  the 
values  of  the  parameter  shown  in  Table  1.  Test  cases  were  integrated  for  an  adiaba- 
tically  collapsing  bubble  and  were  compared  to  other  investigator's  results  to  pro¬ 
vide  an  initial  check  of  the  integration  technique.  To  serve  as  a  basis  for  comparison, 
in  each  case  the  collapse  of  a  vapor  bubble  was  computed  with  no  heat  release  in 
the  vapor  phase  and  shown  as  the  Q=0  curve  in  Figures  1  through  4.  The  curves 
of  pressure  and  temperature  histories  in  the  bubbles  during  collapse  are  plotted  as 
functions  of  the  bubble  radius  which  varies  in  time  very  nearly  like  the  radius  of 
a  collapsing  adiabatic  bubble.  As  shown  in  Figures  1  and  2,  the  temperature  (and 
pressure)  in  the  bubble  during  collapse  follow  the  non-reacting  case  until  a  critical 
temperature  is  reached  where  a  rapid  increase  in  temperature(pressure)or  ignition 
takes  place.  Figure  3  shows  the  effect  of  different  overpressures;  the  rate  of  temp¬ 
erature  increase  at  the  higher  overpressure  is  greater  because  less  heat  is  lost  to 
the  liquid  during  the  more  rapid  collapse.  Ignition  at  both  values  of  overpressure, 
as  would  be  expected,  occurs  at  about  the  same  temperature.  The  effects  of  chang¬ 
ing  initial  bubble  radius  is  shown  in  Figure  4.  The  larger  bubble  ignites  at  an 
earlier  stage  of  collapse,  probably  because  there  is  a  lower  rate  of  heat  transfer 
per  uifit  mass  of  vapor  in  the  larger  bubble  due  to  the  smaller  surface  to  volume  ratio 
of  the  larger  bubble.  Thus,  in  general,  larger  bubbles  and  higher  overpressure  lead 
to  earlier  ignition  (in  terms  of  bubble  radius). 

Once  bubble  ignition  occurs ,  the  integration  of  the  governing  equations 
is  terminated  because  the  properties  change  too  rapidly  for  the  program  to  handle. 

In  interpreting  these  results,  it  should  be  realized  that  several  of  the  assumptions 
made  in  deriving  the  governing  equations  are  not  valid  during  the  latter  stages  of 
the  bubble  collapse.  For  example,  as  the  bubble  wall  reaches  a  very  high  velocity, 
the  assumptions  of  an  incompressible  liquid  phase  and  uniform  bubble  interior  are 
not  valid.  Also,  the  assumptions  of  equilibrium  vaporization  and  reaction  kinetics 
are  no  longer  appropriate  for  high  rates  of  change  of  the  bubble  properties.  These 
results  should,  therefore,  be  viewed  as  a  qualitative  indication  of  the  effects  noted 
above. 
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The  high  pressures  generated  by  the  ignition  process  may  cause  a  reversal 
of  the  bubble  collapse  and  a  re-expansion  of  the  bubble.  The  question  then  arises 
as  to  whether  the  bubble  expansion  would  extinguish  the  reaction  or  whether  the 
reaction  would  continue  to  drive  the  bubble  past  its  original  size.  If  the  latter 
occured,  this  could  be  taken  as  a  criterion  for  the  ignition  of  a  bulk-cavitated 
liquid  containing  many  bubbles. 

The  solution  of  the  rebounding  bubble,  which  requires  that  liquid  com¬ 
pressibility  and  finite  rate  kinetics  be  included  in  the  model  may,  therefore,  lead  to 
an  Ignition  criterion  for  LVD  and  other  hot-spot  ignition  phenomena. 
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Figure  1. 


Figure  2. 


Figure  3 . 


Figure  4. 


Figure  Captions 


Collapse  of  vapor  bubble  with  overpressure  of  107  dynes/cm2  and  in- 
tial  radius  of  0ol  and  showing  effects  of  different  activation  energies 
compared  to  bubble  collapse  with  no  heat  release  or  reaction  (0=0). 

(a)  Temperature  at  the  bubble  wall  vs.  bubble  radius 

(b)  Pressure  in  the  bubble  vs.  bubble  radius 


Same  calculation  as  in  Figure  1  but  with  a  higher  overpressure,  109 
dynes/cm8 . 

(a)  Temperature  at  the  bubble  wall  vs .  bubble  radius 

(b) .  Pressure  in  the  bubble  vs.  bubble  radius. 


(a)  Bubble  wall  temperature  and 

(b)  pressure  in  the  bubble  during  collapse  compared  for  the  two  different 
overpressures. 


(a)  Bubble  wall  temperature  and 

(b)  pressure  in  the  bubble  during  collapse  compared  for  two  different 
initial  radii. 
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Figure  3b,  Pressure  In  The  Bubble  During  Collaspe  Compared  For  The 
Two  Different  Overpressures. 
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109  dynes/cm3 
45 


Figure  4a,  Bubble  Wall  Temperature 


